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ABSTRACT 

In this study I present new K and L'-band infrared photometry and 2-4//m 
spectra of ten Class I and flat spectrum stars forming within the Taurus dark 
cloud complex. Nine sources have H2 ^=0-1 S(l) emission, and some show mul- 
tiple H2 emission features in their K-band spectra. Photospheric absorptions 
characteristic to low mass stars are detected in five of the targets, and these 
stars were fit with models to determine spectral type, infrared accretion excess 
veiling (r^^ and tli) and dust temperatures, estimates of visual extinction and 
characteristics of the 3/im water-ice absorption. On average, the models found 
high extinction values, infrared accretion excess emission with blackbody tem- 
peratures in the 900-1050K range, and 3/um absorption profiles best fit by water 
frozen onto cold grains rather than thermally processed ice. Five techniques 
were used to estimate the extinction toward the stellar photospheres; most gave 
vastly different results. Analysis of emission line ratios suggests that the effect 
of infrared scattered light toward some protostars should not be neglected. For 
stars that exhibited Br7 in emission, accretion luminosities were estimated us- 
ing relations between L^cc and Br7 luminosity. The young stars in this sample 
were preferentially chosen as variables, but they do not have the accretion dom- 
inated luminosities necessary to put them in their main stage of mass-building. 
The characteristics of the 2-4/im spectra are placed in the context of existing 
multi-wavelength data, and five of the stars are catagorized as reddened Class lis 
or stars in transition between Class I and II, rather than protostars embedded 
within massive remnant envelopes. 

Subject headings: stars: formation, stars: pre-main sequence, ISM: dust, extinc- 
tion, infrared: stars 
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1. Introduction 

Ever since the identification of T Tauri stars as young solar analogs (Ambartsumian 
1947; 1949), astronomers have sought to better understand the youngest phases of stellar 
evolution. The youngest optically revealed stellar objects (YSOs) are still obscured by their 
natal envelopes and ambient cloud material, and hence are very faint in visible light. Un- 
derstanding the stellar and circumstellar properties of these young stars has proven feasible 
with the current generation of sensitive instrumentation and large aperture telescopes. The 
Class I protostars represent the youngest YSOs which are presently identified at optical and 
near-infrared wavelengths. 

In the accepted classification scheme for early stellar evolution. Class I, or protostellar 
T Tauris, are embedded young stars that are intermediate between extremely young stellar 
embryos (Class 0) and the more familiar Classical T Tauri Stars (CTTSs, Class lis; Adams, 
Lada & Shu 1987; Lada 1987). The definition of these classes came from the shape of the spec- 
tral energy distribution (SED) from 2-25/im. The spectral index, a. = dlog(AFA)/dlog(A), 
represents the shape of the SED. Sources that have a >0.0 are designated Class I with rising 
SEDs in the infrared, and those with a ~0.0 are fiat spectrum objects (Adams, Lada & Shu 
1987; 1988). The SED shapes of Class Is are well modeled as star-|-rotating disk systems with 
massive infalling envelopes and outfiows (Kenyon et al. 1993a; 1993b, Kenyon & Hartmann 
1995; Whitney et al. 1997; 2003; Terebey et al. 2006). It has become apparent from both the 
models and observations that a degeneracy in the spectral index exists between embedded 
protostars and more evolved but obscured T Tauris, e.g., CTTSs with disks viewed edge-on 
(Whitney et al. 1997; 2003; Chiang & Goldreich 1997; 1999; Hogerheijde & SandeU 2000). 
Investigations in optical through millimeter regimes have sought to better define the differ- 
ences between Class I stars and reddened Class lis (Motte & Andre 2001; Park & Kenyon 
2002; White & Hillenbrand 2004; Andrews & Wilhams 2005; Eisner et al. 2005). 

One of the outstanding discrepancies between the theory and observation of Class I 
protostars is the "luminosity problem" first described by Kenyon et al. (1990). If mass in an 
infalling protostellar envelope is accreting efficiently onto the stellar surface instead of only 
onto the circumstellar disk, then the luminosity from mass accretion should dominate the 
stellar photosphere by a factor of ten or more. However, for the Class I stars that have been 
investigated, the mass accretion luminosities do not differ significantly from those of Class 
lis (GuUbring et al. 1998; MuzeroUe et al. 1998; White & Hillenbrand 2004). This has led 
to speculation that either Class I protostars are not in their main accretion phase and their 
mass was accumulated at a younger epoch (White & Hillenbrand 2004), or that eruptive 
events of enhanced mass accretion in FU Ori-like outbursts are responsible for building the 
bulk of the stellar mass (Kenyon et al. 1990; Calvet et al. 2000). While the fiux of some 
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Class Is is known to vary (Leinert et al. 2001; Park & Kenyon 2002; Forbrich, Preibisch & 
Menten 2006), little is known about the cause and character of variability in most embedded 
protostars. 

Spectral characteristics in the near-infrared provide important opportunities to study 
the environments and natures of YSOs (Greene & Lada 1996). The K-band (2.0-2.4/[/m) 
spectra of Class Is typically have continuum emission that is flat or rising steeply toward 
longer wavelengths, and are often featureless aside from emission in atomic or molecular 
hydrogen transitions which indicate ongoing mass accretion and outflow activity, respectively 
(Greene & Lada 1996; 1997; Prato & Simon 1997; Muzcrollc ct al. 1998; Davis et al. 1994; 
Gomez et al. 1997). On average. Class Is have greater K-band infrared accretion excesses 
than the more evolved CTTSs, and stellar photospheric absorption features characteristic 
to their spectral types are often veiled to undetectability by the excess emission (Greene & 
Lada 1996; 1997; Luhman & Rieke 1999; Doppmann et al. 2005; Ishii et al. 2004). 

2-4//m spectroscopy provides many useful indicators to determine spectral type, accre- 
tion and outflow activity, and the characteristics of material along the line of sight to young 
stars (Whittet et al. 1988; Greene & Lada 1996; 1997; Prato & Simon 1997; Teixeira & 
Emerson 1999; MuzeroUe et al. 1998; MuzeroUe et al. 2003). Yet, this wavelength region 
has not been used extensively for Class Is. In this work, I present 2-4/xm spectroscopy of a 
sample of ten protostars in the Taurus dark cloud complex. These stars were chosen because 
they had Class I or flat SED shapes (Kenyon & Hartmann 1995) and had evidence for near 
infrared flux variability (Park & Kenyon 2002). In §2 the observations arc discussed. In §3 
I present the spectra and the models used to characterize the stars and §4 and §5 discuss 
extinction and accretion properties, respectively. In §6 I place the results in the context of 
existing data and briefly discuss each star. 

2. Observations and Data Analysis 

The data for this project were obtained at the NASA Infrared Telescope Facility (IRTF) 
3 meter telescope at Mauna Kea on UT 2003 Dec. 5 and 6. K and L'-band images and 2-4/im 
cross-dispersed spectra were acquired. During the observations, the sky was photometric and 
the seeing was 0."6-0."8. The facihty near infrared imager and spectrometer "SpeX" was 
used to take the data (Rayner et al. 2003). SpeX consists of a 1024x1024 InSb array, 
Bigdog, which has a plate scale of 0."15/pixel and is used for spectroscopic observations, 
and a 512x512 InSb array, Guidedog, which has a 0."12/pixel plate scale and serves as a 
guiding/imaging camera. The targets for this project were chosen because they are potential 
variable stars; they showed >0.15 mag of deviation from their mean K-band magnitudes in 
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the study of Park & Kenyon (2002). Most of the young stars were not spatially extended by 
more than FWHM~2."2 (2 pixels) in the study of Park & Kenyon (2002). 

Table 1 presents the list of the stars observed for this project. Also included in Table 
1 are the RA and Dec. for the targets, the SED classifications, the standard deviation 
of multiple measurements of the magnitude, ax, which provides an estimate of the stellar 

variability level, and the K and L' magnitudes determined from the data presented in this 
work. (Kenyon & Hartmann 1995; Park & Kenyon 2002; White & Hillenbrand 2004). For 
this study, stars with 2-25/im spectral indices of less than 0.3 are identified as Flat spectrum 
sources rather than true Class Is (Table 1; Adams, Lada & Shu 1988). Photometry was 
obtained for IRAS 04108+2803A because it hes in the same image as 04108+2803B, but no 
spectroscopy was acquired for this source because its value for ax did not classify it as a 
potential variable star (Park & Kenyon 2002). 

2.1. Imaging 

For the 11 stars in Table 1, images were obtained using the Guidedog camera with the 
K and L' filters. The integration times were 160 seconds at K and 80 seconds at L'. This 
total time was derived from summing 40 images, each individual image was an 0.8 second 
integration with 5 coadds in the K-band, and 0.2 second integration with 10 coadds at L'. 
For sky subtraction, the data were obtained by cycling through an a/b beam-switch pattern 
using telescope offsets of 15". Photometry for each of the sources was done using a 5" circular 
aperture. The photometric standard 16 Tau was observed 3 times per night for K-band fiux 
calibration and 2 times for the L'-band; multiple photometric standards were not observed. 
Estimated uncertainties arc ~0.08 and ~0.2 magnitudes in the K and L'-bands, respectively. 
Columns 6 and 7 of Table 1 present the final K and L' magnitudes of the stars. 

2.2. Spectroscopy 

The spectra were obtained using the Bigdog array with a 0."5 slit in the long wave- 
length cross-dispersed setting, LXD1.9, which provides simultaneous spectral coverage of 
the 1.9-4. 1/im region. In order to correct atmospheric features at 3.0-3.4/im, the spectra 
were obtained by observing multiple 30 second exposures with a beam-switch of 7."5 that 
kept the star in the slit at all times. For the success of telluric correction, observations of an 
AO spectral type star were made within <20 minutes in time and less than 0.2 difference in 
airmass of each science target. 
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The data were reduced using the "SpexTool" data reduction package provided by the 
IRTF (Gushing et al. 2004). The spectra were extracted using an optimal extraction with 
a PSF fitting radius of 1" and an extraction aperture of l."2. The SpexTool tasks xtellcor, 
xmergeorders and xcleanspec were used for telluric correction, to merge the orders into a 
single spectrum, and to clean bad pixels. The photosperic hydrogen features present in 
the AO spectral telluric calibrator were removed using the "xtellcor" data reduction routine 
that uses the deconvolution method described by Vacca et al. (2003). In all cases, good 
correction of the hydrogen photospheric features in the AO stars was obtained using this 
method. To test the success of the telluric correction, multiple AO stars were run through 
"xtellcor" as the science targets. Correction of the spectra for these stars was done using 
AO stars observed nearby in time and airmass. The resulting telluric correction was good 
in the 2.40-2.44/im region where H2 Q-branch emission in the spectra of young stars can be 
adversely affected; remnant telluric features were typically not seen above the noise in the 
continuum in the AO stars. The fluxes derived from the K and L' images of the young stars 
were used to calibrate the spectra. 

Observations of thirteen late type dwarf and giant stars with spectral types of Kl 
through M6 were acquired for comparison with the spectra of the young stars. These spectral 
type templates were necessary for use in spectral model fits. The stars were chosen based on 
their spectral types listed in the Hipparchos catalog, and their brightness in the infrared (all 
had K<7.5). These spectra were acquired either on UT 2003 Dec. 5 and 6 at the same time 
as the observations of the young stars, or in a previous session of observing for a separate 
project (UT 2002 May 30 - June 1). Observations of the spectral templates were acquired 
with the same instrument configuration as the young stars, with three to five a/b beam- 
switched cycles (depending on the brightness of the star) of 30 second individual exposures. 
The spectra were extracted and corrected for telluric features following the method described 
above for the young stars. Additional spectral type comparison stars were also obtained from 
the SpeX library of stellar standards available online through the IRTF webpage (Gushing 
et al. 2003)0. 

The final spectra of the Class I and fiat spectrum stars have resolution R ^ 1100 in the 
K-band. Figure la through Ij present the 1.95-4.07/Ltm spectra of the 10 protostars and the 
lower panels show a zoomed view of the K-band region. 



^Spectral template standards obtained with the SpeX instrument are available at: 
http://irtfweb.hawaii.edu/spex 
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3. Results 

For most of the targets, the measured K-band magnitudes are comparable to those 
found from past studies and vahies reported by 2MASS. A few stars are 0.1-0.2 magnitudes 
brighter than 2MASS, but this is perhaps expected from the profile of the shorter wavelength 
Kg filter used to acquire the 2MASS data and the reddened spectral shape of most of the 
targets. One source, IRAS 04239+2436, has brightened to 9.88 mag compared to the average 
of 10.33 from Park & Kenyon (2002). This was one of the more variable stars in the sample 
(Table 1), but given the inherent uncertainties the photometry presented here and that of 
Park & Kenyon (2002), the variation in magnitude for IRAS 04239+2436 is only at a ~3cr 
level of significance. It is also difficult to reach any new conclusions on the variability of 
these stars because of differences in filter profiles, photometric apertures, and plate scales 
used for our respective works. None of the stars were spatially resolved with respect to the 
images of the flux standard star. The spatial pixels were too small and the total integration 
times too short to examine with confldence the character of low surface brightness nebulosity 
in the vicinity of these young stars. 

In Figure la through Ij, the spectra of the stars are plotted; the upper panel presents 
the 2-4//m region and the lower shows a close-up of the K-band. Table 2 provides a summary 
of the emission and absorption features in the 2-4yum region that have been detected in these 
spectra and are used to better understand young stars. Also included in Table 2 is a note on 
the information the line species provides. Young, unveiled, low mass stars have absorption 
features in their K-band spectra characteristic to their spectral type (Wallace & Hinkle 1997; 
Greene & Lada 1996; e.g. Figure lb). Photospheric absorption features are detected in five 
of these stars. 

In low mass protostars, molecular hydrogen emission typically arises from regions where 
outflows plow into the ambient ISM and shock the local material. Emission from H2 v=l-0 
S(l) at 2.12yum was detected in all but one of the stars (IRAS 04181+2654B), and a few 
stars showed emission from multiple molecular hydrogen transitions in the K-band (Figure 
Ig). Unfortunately, the signal-to- noise (S/N) of the L-band spectra were too low to detect 
H2 features present in this range. Additionally, the H2 v^2-l S(l) emission at 2.24y[/m was 
not seen in any of the spectra. Atomic hydrogen emission arises from the magnetospheric 
accretion of material onto the central star, features in the infrared are surrogates to the 
more familiar Ha emission in the optical (e.g. Figure le; Prato et al. 1997; MuzeroUe et 
al. 1998). In addition to Brj in the K-band spectra, strong emission from Bra (4.05//m) is 
also detected in many of the stars. Table 3 presents the line fluxes of the emission features 
detected in the spectra and 3cr detection limits are included where the emission features were 
not detected. 
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Beyond a certain threshold extinction, ices frozen onto interstellar or circumstellar dust 
grains give rise to the solid-state absorptions from 0-H and C-H stretch-mode oscillations 
at 3.05 and 3.47/xm, respectively (e.g. Figure li; Whittet et al. 1988, Teixeira & Emerson 
1999, Brooke et al. 1999). The 3yum water ice absorption in the L-band spectra are detected 
at varying levels of confidence in all but one of the young stars (IRAS 04264+2433). This 
absorption can arise from ices frozen onto grains in the circumstellar disk, natal envelope, 
or from ambient material in the Taurus dark clouds. The shape of the water ice absorption 
profile is temperature dependent, a narrow, more sharply peaked shape is indicative of warm, 
annealed ice and a broad, rounded profile traces cooler dust (Gerakines et al. 1995; 1996). 
Investigating the profile shape of the water feature can provide information on thermal pro- 
cessing in the physical environment where the absorption occurs. Whittet et al. (1988) 
showed that a relation exists between water-ice absorption optical depth (rjce) and visual 
extinction toward young stars embedded within and observed through the Taurus star for- 
mation region. Further work by Teixeira & Emerson (1999) showed that the column density 
of water is a better tracer of the visual extinction toward obscured T Tauri and background 
stars. 

Table 4 presents parameters that characterize the water-ice absorption seen in the spec- 
tra. Column 1 is the blackbody temperature that best fits the spectral continuum shape in 
the 2.0-2.3 and 3.5-4.0//m regions. The additional values described in Table 4 were derived 
by dividing by the blackbody and removing the continuum shape from the spectra. For IRAS 
04181-I-2654B no blackbody temperature was found to fit the K and L-band continuum with 
less than 15% average deviation, so a 2nd order polynomial fit was used. The second and 
third columns in Table 4 are the wavelength of the peak absorption in the ice feature, Ajce, 
and the FWHM of the water ice absorption, Auice- Columns 4 and 5 are the optical depth, 
Tice, and the ice-band column density, Nj^e determined from the data, respectively. Water ice 
column density is found using the relation Nj^e = Tice^^/-^ where t:,,,,, is the optical depth of 
the water ice, Au is the FWHM of the feature in cm^^ and A is the integrated absorption 
cross section, A= 2.0 x 10"^^ molecules/cm^ (d'Hendecourt & AUamandola 1986). 

Uncertainties in the the FWHM of the water ice absorption can be large because of 
difficulty determining the water ice profile shape in the region of telluric absorption (2.55- 
2.9//m). Young stars often have warm material in the inner regions of their circumstellar disks 
where water-ice has vaporized, so the 3/um absorption does not trace all material along the 
line of sight to the stars. Because of these factors, the uncertainties in the visual extinction 
derived using the relations from Teixeira & Emerson (1999a) or Whittet et al. (1988) can 
also be large. As a result, relations between water-ice optical depth or column density and 
the visual extinction are not well calibrated for Class I protostars. 
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3.1. 2-4//m Spectral Synthesis Models 

To better determine the natures of the stars, a 2-step method for spectral synthesis 
modehng was used. First, models were fit to the K-band spectra to determine the stellar 
type and infrared veiling at 2.22/im. The resulting information was then used to generate 
a more complete model for visual extinction, ice absorption and infrared excess over the 
full 2-4//m range. The underlying spectral type of the star is crucial for determining an 
accurate model, so fits were done only for the five stars that had photospheric absorption 
features characteristic to low mass stars. IRAS 04264+2433 has modest infrared veiling and 
an estimated spectral type of Ml±2 (Doppmann et al. 2005; White & Hillenbrand 2004), 
but line ratio analysis (§4) reveals that this star has a strong scattered light component even 
at infrared wavelengths. I do not see infrared photospheric absorption features in this star, 
and since strong scattering would underestimate the extinction to the stellar photosphere I 
do not fit spectral models to it. Models were also not generated for the other four stars that 
did not have spectroscopically determined stellar types. 

3.1.1. The K-band Model: Spectral Type and 2.22iim Veiling 

The K-band spectra for the 5 young stars with photospheric absorption features were 
used to determine the stellar spectral type, the infrared accretion excess (veiling) at 2.22//m, 
and an initial estimate of the visual extinction. To model these spectra, I used a min.- 
imization code written in the IDL programming language which was kindly shared by L. 
Prato and is described in detail in Prato, Greene & Simon (2003). The routine was modified 
for use with the data obtained for this study. The spectral type templates used for the 
analysis were from observations that I made for this purpose, as well as data from Gushing 
et al. (2003) that are available onhne. 

The young star spectrum was fit to a series of models of a spectral type reference star 
which had been veiled with an infrared emission excess and then reddened by extinction 
using an ISM extinction law (A(A) = A^[0.55/A]^-^; Prato , Greene & Simon 2003). The 
K-band veiling, Yk, is defined as the ratio of the infrared excess at 2.22yum to the intrinsic 
stellar continuum flux: F/^(2.22y[xm)/Fj„((2.22yum). The excess is assumed to be wavelength 
independent. A search for the infrared veiling value and a best-fit visual extinction was 
done for a given spectral type standard. Typical models searched over values of r^^O.O to 
1.5, and through a range of 20.0 mag in A^ (typically from A^=10.0 to 30.0). The search 
was repeated for several different stellar spectral type templates, and the final model was 
adopted from the star which gave the best fit in the least squares sense to the continuum 
shape and photospheric absorption features of the YSO. The wavelengths where H2 and Br7 
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emission is pronounced were not used in the range for the fit. 

The spectral type and infrared veihng values at 2.22//m that best-fit the data are pre- 
sented in the first two columns of Table 5. Figure 2 shows the K-band spectra for each star, 
overplotted in blue is the best veiled and obscured spectral model. The models give good to 
excellent fits to the continuum shape and photospheric absorption features of Mg (2.28//m), 
Na (2.20yum), Ca (2.26^,m) and the CO band-heads (2.29-2.38Atm), even though the latter 
features can also be affected by non-photospheric absorption or emission in young stars. In 
all cases, dwarf spectral type standards gave better fits than giant stars to the continuum 
and absorption features. For IRAS 04181+2654A, 04181+2654B and 04295+2251 this anal- 
ysis is the first spectroscopic determination of a stellar type for the underlying photospheres 
(though Doppmann et al. 2005 estimate a temperature of 3400K for 04295-1-2251). For 
the other two stars, the spectral type determined here is consistent with previous estimates 
(White & Hillenbrand 2004). 



3.1.2. The 2-4iim Spectral Model 

Once the best K-band models were generated, the spectral type and tks were used 
to generate more complete models over the full 2-4yum range. The search over the full 
wavelength region allowed for a more accurate determination of the extinction and a fit of 
wavelength dependence to the excess emission. MuzeroUe et al. (2003) showed that infrared 
excess emission in the 2-5/um spectra of Classical T Tauri stars is consistent with blackbody 
emission with a single characteristic temperature. The available laboratory data of water- 
ice absorbance also allows for fits of the 3/im profile shape to search for variations in ice 
temperatures in the circumstellar environment. 

The synthesis model constructed for this test was adapted from the K-band spectral 
code to use a wavelength dependent infrared excess and an optical depth model which had 
ice profiles included. The observed 2-4yLim spectrum for each protostar is defined as: 

Fobs^ C[(F,„t+k(A))e--W] 

Where F^nt is the intrinsic photospheric fiux of the young star estimated from the spec- 
tral standard. k(A) is the wavelength dependent accretion excess, which is described using 
blackbody radiation with a single temperature, Tfi^st, constrained by the required value 
at 2.22;um determined by the K-band models (MuzeroUe et al. 2003). The wavelength 
dependent optical depth, r(A), has two components: one from the ISM extinction law as de- 
scribed for the K-band model in §3.1.1 (Prato, Greene & Simon 2003), and a second optical 
depth component, Tjce, was derived directly from the laboratory absorbance data for water- 
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ice (Gerakines et al 1995; 199(u). The absorption profiles seen in the stars were compared 
to water ice at lOK, 50K, 80K and 120K to search for thermal processing in the absorbing 
material The parameter C is a scaling constant to normalize the models to the flux for each 
star at 2.22/xm in the calibrated spectra (Figure 1). A minimization search was done over 
the three parameters, A^,, Tdust and Tjce for each of the stars in Table 5. 

Early in the modeling process it became apparent that all of the targets exhibited the 
known absorption excess in the long wavelength wing of the 3/xm feature (Van de Bult et 
al. 1985; Smith et al. 1989). To date, it is still unclear what process or species causes this 
excess absorption. Explanations for its existence could include nitrogen or carbon containing 
molecules frozen onto the water-ice mantles, affects of the profile shape from non-spherical 
grains and grain orientation, or overlaid absorption of a broad water vapor feature in the 3/xm 
region (Smith et al. 1989; Zinov'eva 2005). Regardless of the cause of the excess absorption, 
the success of the model fits required the 3.25 to 3.75/im wavelength region to be excluded 
in the analysis (consistent with the results of Smith et al. 1989). The wavelengths with 
atomic and molecular hydrogen emission were also excluded in the fit. The final models 
provided the best simultaneous fit to the young star spectra in the ~2.0-2.4/im, 2.9-3.25/im 
and 3.75-4.04|Um wavelength regions. 



3.1.3. Model Results 

The last four columns of Table 5 present values for A„, r^/, Tdust and r^ce from the best 
model for each target, yl' is the L-band accretion excess at 3.76;um, defined as ri/=F/R(3.76/im) /Fj„t(3.76/ii] 
and is a natural result of this spectral fit. For all stars, L-band accretion excess values were 
3.5 to 5 times greater than those derived for the K-band. Figure 3 shows the 2-4/im spectra 
of the targets with the best model overplotted in red. The visual extinctions from the best 
model were within ~2.5 magnitudes of the values determined from the K-band model. In all 
but one case, the models give a good fit to the K-band spectra, the 2.9-3.25//m wavelength 
region, and wavelengths greater than 3.8;Ltm. 

The uncertainty in a given model was dependent predominantly upon A^ and Tdust- 
Changes in the continuum model shape did not strongly affect the fits for Tice', variations of 
only ~ ±0.07 in the optical depth were typical for a range of extinctions and dust temper- 
atures. For different normalization constants, several values for A^, and T^^^^ could fit the 
continuum shape (Aspin 2003). However, uncertainties in A„ and Tdust can be quantified by 
investigating the reduced surface for a single value of C which normalized the best model 



^Laboratory water ice absorbance data are available online at: 



http: / /www. strw.leidenuniv.nl/ ~lab 
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to the flux calibrated young star spectrum (Prato, Greene & Simon 2003). Figure 4 presents 
the contour plot for the models to IRAS 04295+2251. The contours of x^=l-01, 2, 3, 6 
and 10 are included. Overplotted in blue in Figure 3e for IRAS 04295+2251 is a model that 
has a reduced value of 2. The la uncertainties for the 2 parameter models are defined 
by the x^=3.2 surface (approximated by x^=3; Prato, Greene & Simon 2003). Figure 4 
shows that the contours of the surfaces are narrower for A^ than Tdust- The models in 
the 2-4/xm region were more dependent upon changes in extinction than dust temperature. 
In most cases a narrow range of 2-3 A^ was implied by the best fits. 

No model for IRAS 04181+2654B could be found which fit the K and L-band continuum 
shape simultaneously. The "best" model presented in Figure 3c and plotted in red deviates 
from the spectral energy in the K-band continuum by ~15%, on average. Also overplotted in 
the figure in green is the x^ fit only to the K-band and 2.9-3.25/xm region. All models which 
fit the steep slope to the K-band continuum overestimated the longer wavelength fiux by a 
factor of two. The value for ik was 0.0 from the K-band model of IRAS 04181+2654B, so 
even in the absence of infrared veiling emission this object can not be well fit from attenuation 
by ISM extinction alone. Discussion of the spectra and the nature of this strange source is 
continued in §6.1. 

The blackbody dust temperatures that fit the wavelength dependence of the infrared 
excess are in the range of ~900-1050K. This is warmer than the blackbody fits to the con- 
tinuum in the raw data (Table 4), but not as high as the 1100-1400K temperatures derived 
for infrared excesses in the Glassical T Tauris studied by Muzerolle et al. (2003). Although 
infrared dust temperatures of greater than 1200K could not be excluded by the models, 
excesses with temperatures less than ~650 were always ruled out with ~3cr confidence. The 
visual extinctions determined for the spectral models are typically greater than values deter- 
mined using other techniques (§4). The continuum shapes and deep profiles for the ice-band 
absorption are relatively well-fit by the models for IRAS 04181+2654A and 04295+2251. 
However, the shallower absorptions seen in IRAS 04158+2805 and 04260+2642 are fit more 
poorly (see further discussion in §6.1 for IRAS 04158+2805). The ice absorption profiles 
were always better fit by optical depth models for cold ice (Tj^e < 50K) versus warmer, 
annealed ice (>80K), indicating that significant thermal processing of the icy material along 
the fine of sight toward these cool dwarfs has not occurred. 

The 2-4/im modeling method provides a powerful way to determine the properties of 
young stars when the spectra can be well described by a veiled and reddened photosphere. 
Yet, a significant fraction of the infrared fiux from Glass I protostars is scattered off the 
surface of the disk or envelope (Kenyon et al. 1993b; Whitney et al. 1997; Park & Kenyon 
2002). The infrared images with high spatial sampling can still be dominated by scattered 
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light (Padgett et al. 1999; Terebey et al. 2006). Sellgren et al. (1996) showed that in 
the absence of PAH emission in the 3.3^m region, infrared reflection nebulosities can be 
approximated using Rayleigh scattered light with a wavelength dependence. Ishii et al. 
(2004) modeled their spectra of an infrared reflection nebulosity near an embedded protostar 
with some success. As a test, a very simple scattered light model with an isotropic phase 
function and Rayleigh wavelength dependence was introduced into the 2-4//m spectral fits. 
This was accomplished by attenuating the standard plus infrared excess model, Fj„t+k(A), 
by Rayleigh scattered light before the star was reddened with ISM extinction. This rudimen- 
tary investigation showed that if the contribution of scattered hght is appreciable (~50% of 
the flux at 2.2//m), the spectral models that do not take scattering into account will under- 
estimate both the visual extinction and the infrared veiling. As a result, the values for A^ 
and tk determined from this analysis are probably lower limits for the zero scattered light 
case. 



4. Extinction Toward Embedded YSOs 

Accurate values for the visual extinctions are needed to derredden the observed fluxes 
and derive underlying stellar luminosities, but determining the obscuration toward an em- 
bedded protostar is challenging because of the unknown contribution from scattered light. 
The available infrared photometry coupled with the 2-4/im spectral data provides five inde- 
pendent opportunities to measure visual extinctions toward these young stars (Table 6): 

1) A„(ice) - Determined by equating the column density of water ice along the line of 
sight to an A^, using the relation presented in Teixeira & Emerson (1999). The applicability 
of the Nice to At, relation questionable; the scatter is large for embedded sources with A^ >20 
determined from infrared photometry. The column density will not trace warm circumstellar 
material where ices have sublimated off of dust grains (T>210K). Hence, the extinction 
derived from water ice absorption is a very rough estimate of the amount of obscuring 
material along the line of sight to the embedded photosphere. 

2) A^(phot) - Estimated by dereddening J, H, and K-band magnitudes from 2MASS to 
the Classical T Tauri star locus (Meyer et al. 1997). 

3) At, (spec) - Derived by using the 2-4/xm spectral synthesis models to simultaneously 
determine infrared accretion excess and A^, (see §3.1). 

4) A„(H2) - Found using molecular hydrogen line emission from the v=l-0 Q(3) and S(l) 
transitions which arise from the same state and have an intrinsic ratio of Q(3)/S(l)=0.7 
(Turner 1977). A„ values are determined by assuming an ISM extinction law (defined in 
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§3.1.1) and calculating the attenuation required to redden the 0.7 ratio to the observed value. 
This line analysis provides a robust tracer of material along the line of sight to the emitting 
region, but the molecular hydrogen may not sample material toward the photosphere because 
it could arise from shocks that are spatially extended from the star. 

5) A„(BrQ;/Br7) - Estimated by using Bra to Br7 emission hue ratios. In a Case B 
treatment where Lya photons are optically thick (Baker & Menzel 1938), and for a wide range 
of temperatures and densities (assumed to be applicable to the atomic hydrogen emitting 
environments of young stellar magnetospheres). Bra is intrinsically 2.85 to 3.17 times brighter 
than Br7 (Brockelhurst 1971; Hummer & Storey 1987). A^s arc estimated by assuming an 
intrinsic ratio of ~3 and reddening this to the observed value using an ISM extinction law 
(§3.1.1). 

Table 6 summarizes the visual extinctions for the ten YSOs that were found using these 
different methods. Each relies on assumptions that scattered light does not affect the visual 
extinction, and that the obscuration along the hue of sight follows an ISM extinction law (as 
defined in §3.1.1). Four of the stars have Q(3) and S(l) emission for H2 hue ratio analysis, 
five have spectral model fits, and seven have A^s from Brackett ratios. The independent 
extinction indicators give different estimates for A^. IRAS 04295+2251 is the only star 
with extinctions determined using all methods, and the resulting values differ by more than 
20 mag. Dereddening the infrared photometry is the most widely used technique in the 
literature. 

Extrinsic to the stellar photosphere and circumstellar emitting region around a YSO, 

the infrared flux is attenuated by both extinction and scattering. Emission line ratios can 
provide information on the decrease of flux from scattering compared to reddening. Plotted 
in Figures 5 and 6 arc the H2 Q(3) versus S(l) line fluxes and the Bra versus Br7, respectively. 
In Figure 5, a line showing the intrinsic Q(3)/S(l) value is plotted, as are lines that show 
the behavior of the ratio when the emission is viewed through At,=10mag and 20mag. For 
the Brackett emission (Figure 6), similar lines are overplotted but a range of values instead 
of a single line is identified to take into account the ratios for differing physical environments 
of the atomic hydrogen emitting region (Brockelhurst 1971; Hummer & Storey 1987). For 
intrinsic hydrogen line emissions seen in pure (Raylcigh) scattered light, the H2 Q(3)/S(l) 
ratio becomes 0.41, and the average Bra/Br7 ratio is ~0.24. Lines with these slopes are also 
overplotted in Figures 5 and 6. 

Data points to the upper left of the intrinsic line ratios in Figures 5 and 6 correspond to 

stars that have extinction dominating the scattered light for the ratio, and objects that lie to 
the lower right of the intrinsic line ratio arc dominated by scattered light. A combination of 
reddened and scattered light will make a star appear less obscured in Figures 5 and 6. Most 
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of the sources lie within the range of 0-15 magnitudes in the figures, even though higher 

extinctions of >20 mag are assumed to obscure Class Is and describe their SED shapes 
(Whitney et al. 1997). IRAS 04108+2803B is the only star that has a visual extinction 
determined from the atomic hydrogen line ratio that is greater than for other methods. 

IRAS 04264+2433 has a molecular hydrogen fine ratio of ~0.4, which is less than the 
intrinsic ratio and is consistent with emission viewed in scattered light (Figure 5). This 
source has no detectable Bra emission, but its Br7 flux and Bra detection limit also put it 
underneath the intrinsic line ratio range and in the scattering dominated region of Figure 6. 
This suggests that this star is viewed entirely in scattered light, which is consistent with the 
edge-on geometry adopted to explain observations at optical and mid-infrared wavelengths 
(see also §6.1; White & Hillenbrand 2004; Kessler-Silacci et al. 2005). 

The derivation of multiple A^ values (Table 6) was done to test the applicability of 
techniques used for Class II stars to the more obscured Class I and flat spectrum objects. 
None of the extinctions presented in Table 6 can be deemed correct with a high level of 
confidence because of systematic affects such as infrared excess from disk emission (e.g., 
J, H, K photometry) and unknown levels of scattered hght (Figure 5, 6). Several authors 
are fitting the SEDs of young stars in conjunction with the scattered light images and 
polarization fractions (Whitney et al. 1997; 2003a; 2003b; Eisner et al. 2005; Terebey 
et al. 2006). This method for determining a valid A,, value is promising. The models of 
Whitney et al. (1997) predict that the attenuation of photospheric flux toward an embedded 
protostar can correspond to 20-60 or more magnitudes of optical obscuration. The A^, values 
derived here are at the lower end of this range. In the future, detailed modeling of the multi- 
wavelength SEDs, coupled with polarimetric data and scattered light images, may make 
accurate visual extinctions available for the stars studied here (i.e., Whitney et al. 1997; 
2003a; 2003b; Eisner et al. 2005). 



5. Accretion Luminosity and Mass Accretion Rates 

MuzeroUe et al. (1998) showed that the dereddened Br7 fluxes from Classical T Tauri 
stars correlate with accretion luminosities determined from the hot flux excesses measured 
in the ultraviolet. Accurately dereddening Br7 fiuxes toward these protostars allows for the 
estimation of their mass accretion luminosities. For the five stars that have spectral model 
fits, the Br7 emission is dereddened using the visual extinctions from the best model. IRAS 
04108-I-2803B, 04361-1-2547 and 04365-1-2535 have some of the deepest and most pronounced 
water-ice absorption features in this sample, which suggests strong attenuation of the photo- 
spheric flux. The visual extinctions adopted for these stars are the highest values from Table 



- 15 - 



6. For the remaining two targets, the Br7 fluxes are dereddened using the visual extinction 
from the infrared photometry. 

The visual extinctions used for dereddening the Br7 fluxes are hsted in the flrst column 
of Table 7. Table 7 also presents the dereddened Br7 line flux values, the derived accretion 
luminosities using MuzeroUe et al.'s (1998) relation, the bolometric luminosity (from Chen 

et al. 1995) and the ratio of the latter two. Uncertainties in the accretion luminosities are 
dominated by the ambiguities in the dereddening process (§4). For example, if the Ay is 
underestimated by 5 or 10 magnitudes, which seems plausible based on the values presented 
in Table 6, then the accretion luminosity will be correspondingly underestimated by factors 
of 1.9 and 3.7, respectively. Also included in the lower section of Table 7 are the accretion 
luminosities from the data discussed by MuzeroUe et al (1998) and a comparison to the 
bolometric luminosity. IRAS 04361+2547 and IRAS 04239+2436 overlap in our studies, 
and the accretion luminosities that I derive are 2.5 to 3 times greater than the values they 
present. The sources were chosen for this study because they showed variability in infrared 
flux. Changes of a factor of 2.4 in the observed Br7 flux will cause the accretion luminosity 
to increase threefold (assuming extinction is a constant). Flux variations of this level are 
not surprising for variable YSOs. 

IRAS 04239+2436 is the only star that has an Lacc/LsoZ of >0.5, and three stars have 
Lacc/LsoZ in the 0.2-0.3 range. All but two of the stars are actively accreting material from 
their circumstellar environments, but with the exception of IRAS 04239+2436 the accretion 
luminosities are not a dominant fraction of the total LboI- The stars observed by MuzeroUe 
et al. (1998) that have strong bolometric luminosities typically have lower Lacc/^'Boi fractions 
than the variable YSOs in this sample. 

To relate the accretion luminosity to the mass accretion rate, information on the un- 
derlying stellar mass is required. The discussion in §4.1 of Doppmann et al. (2005) convinc- 
ingly outlines problems using infrared spectroscopy, photometry, and theoretical evolutionary 
models to place embedded protostars directly on a HR diagram to determine their mass. The 
applicability of evolutionary models toward very young stars is uncertain. Moreover, the de- 
rived visual extinctions are known with so little confldence that dereddening infrared fluxes 
to flnd stellar luminosities is also questionable. Nevertheless, models predict that a young 
star evolves vertically in luminosity along the Hayashi track in the 0.5-1.0 million year old 
(Myo) age range (Siess et al. 2000). Thus, I follow the analysis of White & Hillenbrand 
(2004) and estimate a Class I stellar mass based on the location of the one million year 
isochrone at the position of the temperature of the star. The temperatures scale of White & 
Hillenbrand (2004) for stellar spectral types is adopted. The results are shown in the second 
column in Table 8. The two stars that overlap with the White & Hillenbrand (2004) sample 
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have the same mass estimate in this study (IRAS 04158+2905 is substellar). 

GuUbring et al. (1998) estimated that the hberated energy from the accretion is related 
to the mass accretion rate, M, by using a Virial Theorom treatment: Lace— GM*M/R^ x (1/R*- 
l/Kin), where Rj„ is the disk inner radius and has a typical value of 3-5R*. Using the mass 
estimates in Table 8, the corresponding luminosity for a 1 Myo star from the Siess et al. 
(2000) models to estimate radii, and a Rj^ value of 5R*, I have derived the mass accretion 
rates for the five stars in Table 8. 

The accretion rates presented in Table 8 are similar to the average mass accretion values 
of ~lO~^M0/yr for Class II T Tauri stars (Gullbring et al. 1998). These rates are several 
orders of magnitude too low to build the mass of a sun-like star in a ~fcw xlO'' years and 
are not consistent with envelope collapse models (Kenyon et al. 1990; Kenyon et al. 1993; 
Muzerolle et al. 1998; White & Hillenbrand 2004). It is yet unclear if the unexpectedly low 
accretion rates for Class I stars derived by optical and infrared studies are because the main 
mass building phase in a stars life occurs either earher than the Class I stage or by eruptive 
increases in the mass accretion (e.g. FU Ori or EXor outbursts; Kenyon et al. 1990 Calvet et 
al. 2000; White & Hillenbrand 2004). Additional knowledge of the variability of embedded 
protostars would be useful to investigate this. 

IRAS 04158+2805 and IRAS 04260+2642 have mass accretion values derived by White 
& Hillenbrand (2004), these values are also presented in Table 8. Although we used the same 
underlying assumptions for stellar masses and luminosities, the values for M that I determine 
are more than an order of magnitude different from the ones they present. The methods for 

determining these vahics have proven reliable for Classical T Tauri stars (Muzerolle et al. 
1998), so it is tempting to attribute differences in mass accretion rate to the variable nature 
of these stars. However, for IRAS 04260+2642, the limit for the accretion rate that I derived 
is 2.5 orders of magnitude less than the detected value from White & Hillenbrand (2004). 
Since the peak accretion rate enhancement in FU Ori outbursts is estimated to be 2-4 orders 
of magnitude (Calvet et al. 2000), the decrease seen in 04260+2642 seems excessively large 
for stellar variability over the one year time frame that the data for our respective works 
was obtained. In the absence of further information on accretion variations in embedded 
protostars, it seems more probable that the standard techniques used to determine mass 
accretion in different wavelength regimes in Class II stars do not translate well to studies of 
embedded Class I protostars or YSOs observed strongly in scattered light. 
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6. The Nature of the Targets 

Five of the targets included in this study had prior evidence for outflow activity from 
observations reported in the literature, and the other five were not known to have outfiows 
traced by optical emission species, molecular hydrogen emission or CO diagnostics at submm 
wavelengths (Gomez et al. 1997; Bontemps et al.l996; Moriarty-Schieven et al. 1992; Kenyon 
et al. 1998; White & Hillenbrand 2004). Eight stars show significant H2 emission in their 
spectra, and several exhibit emission in multiple transitions (Figure 1). One of the stars 
that does not show appreciable molecular hydrogen emission, IRAS 04365+2535, is known 
to drive a molecular outfiow (Bontemps et al. 1996) and has H2 v—1-0 S(l) emission detected 
at a ~3cr level of confidence by Ishii et al. 2004. In this sample, nine of the ten stars have 
H2 emission. It is not possible to determine conclusively between an origin of the II2 from 
shocks in an outfiow versus UV or X-ray fiuorescent emission. However, the youth and 
low masses and bolometric luminosities of these stars (low UV flux) weigh more toward an 
outflow origin. 

The existing degeneracy in the spectral index between true Class Is and reddened Class 
lis has led a number of authors to try to distinguish the evolutionary states of young stars 
using other methods. In the near-infrared. Park & Kenyon (2002) used the spatial extent 
of surrounding reflection nebulosity to differentiate Class Is from transition or reddened T 
Tauris. Infrared accretion excess emission is usually strong in embedded protostars, and 
values of ick <0.2 are rarely seen for Class Is (Greene & Lada 1996; Luhman & Rieke 1999; 
Doppmann et al. 2005). Kessler-Siliacci et al. (2005) also proposed an evolutionary sequence 
based on the lOyum profile of the silicate absorption and emission; more embedded protostars 
exhibit strong absorption and the more evolved T Tauris exhibit only silicates in emission. 
Though, the use of the sihcate diagnostic alone is ambiguous because source geometry can 
also cause differences in 10/xm profiles (Chez et al. 1991; Chiang & Goldreich 1999; Ciardi 
et al. 2005). Additionally, Motte & Andre (2001) devised an evolutionary sequence based 
on envelope masses in which the younger protostars have resolved mm emission with Meuv 
greater than ~0.1 M©, and more evolved sources have unresolved emission and envelope 
masses less than this. Scattered light levels were determined by Whitney et al. (1997), and 
the embedded protostars were found to have higher infrared polarization fractions. 

Table 9 presents a summary of the observed properties of these ten young stars. The 
first and second columns are the bolometric temperatures and spectral indices, which are two 

parameters that have been used in the past to define the Class Is (Chen et al 1995; White 
& Hillenbrand 2004). Also included are the Lacc/^Boi fractions and water-ice optical depths 
from this work (columns 4 and 5) and data collected from Kessler-Siliacci et al. (2005) and 
Watson et al. (2004) for the silicate emission (column 6), Motte & Andre (2001) for envelope 



-18- 



masses (column 7), and Whitney et al. (1997) for polarization values (column 8). 

Figure 7 presents the spectra for the 10 stars plotted from bottom to top in order of 
increasing K-L' color. Included to the left of the figure is a note on the order of the spectra 
from Figure 1. Seven YSOs have water-ice optical depths greater than one indicating a high 
level of obscuration along the lines of sight, but no apparent relation between Tjce and other 
characteristics is seen. However, trends between the other properties from Table 9 are found. 
The stars that have 10/um silicates in emission typically have lower or unresolved envelope 
masses, lower 'Lacc/'^Boi fractions, and low polarization fractions. These stars are in transition 
between Class I and Class II or arc reddened Class lis. Conversely, the embedded protostars 
have higher average hacc/^Boi fractions, silicates in absorption (where data is available), 
higher envelope masses and stronger polarization fractions. Although their spectral indices 
and bolometric temperatures are consistent with a Class I designation, IRAS 04108-I-2308B 
and IRAS 04264+2433 show many characteristics consistent with a more evolved nature. 

Overlaid on the right side of Figure 7 is a note regarding the nature of these stars from 
the accumulated data (Table 9). IRAS 04108+2308B is plotted as an embedded star in 
Figure 7 to avoid overlapping spectra, but this classification may be incorrect (see discussion 
below). The results presented here are largely consistent with the classifications of these 
YSOs from the other studies. However, I identify IRAS 04181+2654B and IRAS 04295+2251 
as reddened Class lis because they have unusually low infrared accretion veiling emission 
and several indicators of a more evolved nature (sec also §6.1). Observations show that 
40-50% of the stars identified as Class Is from their spectral indices arc more evolved stars 
seen strongly in scattered or reddened light (Whitney et al. 1997; 2003; Motte & Andre 
2001; Park & Kenyon 2002; White & Hillenbrand 2004). In this work, I find the sample of 
ten variable YSOs to be comprised of five true protostars, three in transition between Class 
I and II, and two reddened Classical T Tauris. 

6.1. Notes on Individual Sources 

IRAS 04108+2803B: Strong atomic hydrogen emission features (Br7, Bra, Pf 7) and 
molecular hydrogen ti=l-0 S(l) emission are detected in IRAS 04108+2803B. No absorp- 
tion from photospheric features seen in young solar type stars is found, presumably because 
the infrared excess emission overwhelms the absorptions. This object shows a deep and pro- 
nounced water-ice absorption feature. It has low envelope mass and unresolved mm emission 
(Motte & Andre 2001), low level of polarized infrared light (Whitney et al. 1997), and 10//m 
silicates detected in emission (Kcssler-Siliacci et al. 2005). Chiang & Goldreich (1999) also 
fit the SED of this star as a CTTS viewed through obscuration in its inclined circumstellar 
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disk. 04108+2308B is plotted as a Class I star in Figure 7 (to avoid overlapping spectra), but 
based on the multi-wavelength data listed in Table 9 this placement is probably erroneous. 
It seems more likely that IRAS 04108+2308B is strongly reddened but in transition between 
the Class I and Class II phase, or is an obscured Class II. 

IRAS 04158+2805: Weak atomic hydrogen emission features and molecular hydrogen 

emission arc seen in IRAS 04158+2805. The spectra show photospheric absorption features 
with a best fit spectral type of M6±l. The infrared veiling value of of 0.25 derived for this 
star is significantly less than the 1^=1 -l determined by Doppmann et al. (2005), suggesting 
variable infrared accretion excess emission. 

The 3/im water absorption feature observed in IRAS 04158+2805 is broad and shallow 
with no distinct rounded peak as seen in most other sources. Particularly, the peak ab- 
sorption is at 3.12yum for IRAS 04158+2805 (Table 4), which is shifted with respect to the 
absorption peak of the other stars. The data were not well fit by the laboratory ice-band 
absorbance data (Figure 3a). With a spectral type of M6, this source is potentially a Class 
I brown dwarf (White & Hillenbrand 2004). Late spectral type stars with cool photospheres 
have long been known to exhibit pronounced photospheric water vapor (Woolf 1964). In- 
deed, the broad and poorly peaked absorption profile at 2.9-3.6/Ltm and excess absorption in 
the 2.3-2.5;um region seen in IRAS 04158+2805 is more characteristic of vapor absorption in 
very late type Mira variables than it is of rounded or peaked ice profiles seen in embedded 
YSOs (Elias 1978; Whittct ct al. 1988; Sato 1990). It seems plausible that a pronounced 
photospheric water vapor feature is confusing the identification of water-ice absorption. High 
spatial resolution images show that this source is probably dominated by scattered light at 
2/um (Duchene et al. in prep), so the relative fraction of scattered versus transmissed light 
may also affect the continuum shape in the 2-4//m spectra. The available data (Table 9) sug- 
gest that IRAS 04158+2805 is a late-type star in transition between Class I and II phases. 
Though, White & Hillenbrand (2004) pointed out that the low and unresolved envelope 
mass (Table 9) may be a consequence of the lower central stellar mass, so additional data is 
required to confirm this interpretation. 

IRAS O4I8I+2654A: This star shows weak photospheric absorption features, and a 
best spectral type of M3^^ is found. The K-band accretion excess of rii-=1.3 is the largest 
veiling value derived from these data. Atomic and molecular hydrogen emission features are 
detected, and this star shows a deep water-ice absorption feature. IRAS fiuxes and envelope 
masses are spatially unresolved with IRAS 04181+2654B. The multi-wavlength data shows 
that this star is a Class I. 

IRAS O4I8I+2654B: IRAS 04181+2654B is unique among the spectra in the sample; 
it is the only star that has neither molecular nor atomic hydrogen emission features indi- 
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cating accretion or outflow activity. The K-band region shows a very steep spectral shape 
with photospheric absorption features and a best-fit spectral type of K7±l. Photospheric 
absorption is often veiled to undetectability by infrared disk emission in young protostars. 
Yet, in contrast to the spectra of 041 81 A, B shows no evidence for infrared veiling emis- 
sion. Additionally, IRAS 04181+2654B is undetected in the optical, the estimated I-band 
magnitude is >24, consistent with the very large visual extinction derived from the infrared 
photometry and the K-band spectral modehng (§3). The extinction derived from the K-band 
spectra corresponds within la the extinction determined from the 2MASS photometry. This 
result indicates that there is not a significant amount of scattered light or infrared accretion 
excess emission, both of which affect the A„ determinations from the photometry of Class I 
protostars. 

IRAS 04181-I-2654B shows none of the near-infrared spectral indicators of youth, which 
argues against a Class I nature. Indeed, its protostellar designation is based on spatially 
unresolved IRAS fluxes with 04181A. IRAS 04181-F2654B is too bright in the infrared to 
simply be a background K dwarf or giant star reddened by 40 magnitudes of extinction. A 
K dwarf observed behind the Taurus dark clouds through ~40 mag of obscuration would 
have an estimated K magnitude of ~15, which is 4 mag fainter than 04181B appears (Table 
1). Moreover, its lO^m spectrum has been modeled to have silicates in emission, which is 
an indicator of youth (Kessler-Silacci et al. 2005). Curiously, the L-band flux is 50% fainter 
than expected from a reddened, unveiled young star photosphere (Figure 3c). The enhanced 
extinction and ~35" separation from 04181-I-2654A lead to the hypothesis that 04181-I-2654B 
could be a young star observed through the extended protostellar envelope of A (Motte & 
Andre 2001). If this is the case, the assumption of ISM obscuration may be incorrect and 
this may be why the 2-4/im spectral models could not flt the continuum shape. Further 
spatially resolved observations, particularly at mid-infrared wavelengths, are necessary to 
determine the true nature of this interesting star. It seems most likely that 04181-I-2654B is 
misclassificd as a Class I and is in fact a reddened Class II. 

IRAS 04239+2436: This is the most actively accreting star in this sample. The derived 
accretion luminosity is greater than ~50% of the total bolometric luminosity. It shows 
emission in H2 v=l-0 S(l) and a rich spectrum of atomic hydrogen emission from the BrS, 
Br7, Bra, Pf7, Pfe, and the Hu 14, 15 and 17 transitions. It also exhibits He I (2.06 /xm) 
emission, the CO A v^2 band-heads in emission (also seen by Greene & Lada 1996), and 
evidence for weak emission (~3(t) in the Na features at 2.206/2.208//m. Presence of the Na 
features in emission imply that this star may be more massive than a typical K-M spectral 
type Class I star. Additionally, He I emission is believed to arise from accretion driven winds 
in young stars and its presence here confirms a strong level of mass accretion (Edwards et al. 
2003). The spectra shows deep 3 /im water-ice absorption and a weak C-H stretch absorption 
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at 3.47y[/m. IRAS 04239+2436 is an actively accreting Class I star. 

IRAS 04260+2642: Marginal detection of atomic hydrogen emission features, yet this 
star has molecular hydrogen emission from six transitions in the 1.95 - 2.45/xm spectral 
region. It also shows photospheric absorption features characteristic to young stars and a 
spectral type of K5±2. IRAS 04260+2642 has a low envelope mass and mm emission which 
is spatially unresolved (Motte & Andre 2001). The optical forbidden emission lines from 
this star are narrow, and an edge-on orientation where velocities from the outflow are in 
the tangential direction was adopted to explain this (White & Hillenbrand 2004). IRAS 
04260+2642 is likely in transition between the Class I and Class II phase. 

IRAS 04264+2433: This star shows features from seven different H2 transitions in the 
1.95-2.45|Um spectra. The molecular hydrogen v=l-0 Q(3)/S(l) line ratio is less than the 
intrinsic value (Figure 6), which I interpret as being caused by strong scattered light. IRAS 
04264+2433 is the only star observed for this study that has a 10/im spectrum with silicates in 
emission only (Kessler-Sihacci et al. 2005). The edge-on disk geometry adopted by previous 
studies (White & Hillenbrand 2004; Kessler-Sihacci et al. 2005) is consistent with the data 
presented here. This star is in transition between Class I/Class II and is seen very strongly 
in scattered light at optical and infrared wavelengths. Prior to July 2005, IRAS 04264+2433 
was listed incorrectly in the SIMBAD online database as Elias 6 (Park & Kenyon 2002, 
Motte & Andre 2001). 

IRAS O4295+225I: This star exhibits deep water-ice absorption, atomic and molecular 
hydrogen emission, and absorption features typical of the photospheres of young low mass 
type stars with a best spectral type of K7+2. IRAS 04295+2251 has the lowest infrared 
excess emission detected in this study, rii-=0.09, and a low veiling value was also determined 
by Doppmann et al. (2005). The small infrared excess coupled with the fact that the lOyum 
spectrum has silicates in emission, the infrared polarization is low, and the star has an 
unresolved and low mass envelope leads to the proposition that this is a strong candidate 
reddened Class II star. Eisner et al. (2005) found that the SED and optical/mm images of 
this star are best fit with a disk+obscuration model, although it could not be discerned if 
excess obscuration lies in an extended envelope (Class I star) or could be from ambient cloud 
material (reddened Class II). Density enhancements of ^^20 mag are possible through the 
Taurus star forming region to redden IRAS 04295+2251 if it was located behind the clouds 
(Whittet et al. 1988; Teixeira & Emerson 1999). Yet, unlike for IRAS 04181+2654B, this 
young star is revealed in the i-band (White & Hillenbrand 2004; Eisner et al. 2005). The 
strong visual extinction of A^=28.0 mag determined for IRAS 04295+2251 from the 2-4//m 
spectral model must be caused in part by local disk or envelope material which attenuates 
the infrared light but scatters the optical flux (Eisner et al. 2005). Based on the discussion 



-22- 



presented here and values listed in Table 9, IRAS 04295+2251 may be a reddened Class II 
star. 

IRAS 04361+2547: This object star several atomic hydrogen emission features and 
has a rich molecular hydrogen emission spectrum from seven transitions at 1.95-2.45//m. It 
also exhibits deep water ice absorption and absorption in the C-H stretch mode oscillation 
at ~3.47^m. IRAS 04361+2547 also has the highest envelope mass and greatest level of 
infrared polarization of any star observed for this study. There is no evidence for photospheric 
absorption features. IRAS 04361+2547 is a true Class I protoster. 

IRAS 04365+2535: Strong atomic hydrogen emission features but marginal molecular 
hydrogen emission (~2(t at 2.12//m) are seen toward this source. It also shows continuum flux 
that rises steeply toward longer wavelengths, a very deep, rounded 3 /xm water ice absorption 
feature, and absorption in the C-H stretch mode oscillation at ^3.47 jim. There is no 
evidence for photospheric absorption features and high envelope mass and strong polarization 
confirm its nature as a true Class I embedded protostar. 

7. Summciry 

1) New photometry and 2-4/xm spectra are presented for ten Class I and flat spectrum 
young stars previously found to have variable infrared flux. Molecular and/or atomic hy- 
drogen emission features are detected in all except for IRAS 04181+2654B, and five stars 
show photospheric absorption features characteristic to low mass stars. Evidence for water 
absorption at 3/xm is seen in all but IRAS 04264+2433. 

2) For the five stars that show photospheric absorptions, a two step process was carried 
out to fit models to the 2-4/im spectra. Spectral types, ice absorption characteristics, visual 
extinctions and blackbody fits to the infrared accretion excess were done. For three of 
the stars, the model presents the first spectral type determination, and the other two are 
consistent with past studies. The stellar spectra were best fit by models with cooler accretion 
excess dust temperatures than Classical T Tauris, visual extinctions that were typically 
higher than found using other methods, and water-ice absorption profiles characteristic to 
cold dust grains. 

3) Visual extinctions determined for the young stars using five techniques gave very 
different results. The commonly used method of dereddening infrared fluxes gave A^, vahics 
in the middle range compared to other techniques, and higher values for the extinctions 
were typically found from the 2-4/xm spectral model flts. All of the derivations of extinction 
assumed that the obscuring material follows an ISM law, and that scattered light is not 
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an important factor. This analysis was done to demonstrate the wide range of values 
obtained by applying techniques used for CTTSs to the more obscured Class I and flat 
spectrum protostars. 

4) Extinctions determined by dereddening atomic and molecular hydrogen line ratios 
resulted in low values (0-15mag) for most stars, and negative values two cases. This suggests 
that scattered light does strongly affect infrared fluxes. The hne ratios for IRAS 04264+2433 
are consistent with an orientation where the infrared flux is seen entirely in scattered light 
(such as a disk viewed edge-on). Determining a true value for the photospheric attenuation 
should take scattering into account, and simultaneous modeling of the SED with scattered 
light images or polarimetry provide the most promising technique (Whitney et al. 1997; 
2003; Eisner et al. 2005; Terebey et al. 2006). 

5) Br7 line fluxes were dereddened using a best estimate for the extinction, and accre- 
tion luminosities (and limits) were determined for all stars. As from past works, the stellar 
accretion luminosities are not as high as predicted from envelope collapse models. A differ- 
ence in M of 2.5 orders of magnitude is found for IRAS 04260-1-2642 from this study and 
results from optical spectroscopy (White & Hillenbrand 2004) . This seems excessive for vari- 
able mass accretion for a one-year timescale (§5), and suggests instead that the techniques 
used to derive M at different wavelengths for Classical T Tauri stars may have systematic 
problems when applied to younger stars. 

6) Characteristics of the 2-4//m spectra of the ten stars are placed in the context of 
existing data from the literature. Trends are seen between the observed multi-wavelength 
properties; stars with high envelope masses and polarization fractions typically have higher 
accretion luminosities and silicates in absorption, whereas stars with low, unresolved enve- 
lope masses and low polarizations have evidence for silicates in emission and lower levels of 
accretion. From the accumulated data and new information presented here, two of the stars 
are best classified as reddened T Tauris, and three are stars transition between Class I and 
II and seen strongly in scattered hght. 

I am grateful to Lisa Prato and Saeid Zoonematkermani for sharing their spectral fit- 
ting code for determining stellar spectral types, veiling values and visual extinction from 
K-band spectra of young stars. I thank Lisa Prato and Michal Simon for providing valuble 
suggestions during the course of this project, and Colin Aspin for reading the manuscript 
and offering advice for its improvement. I also thank the anonymous referee who provided 
comments that improved the paper. I appreciate the support provided by the IRTF staff - 
Bill Golish and Bobby Bus for their help with observing, and John Rayner and Mike Gushing 
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Table 1. Target List and New Photometry 



Star Name (IRAS #) 


RA J (2000) 


Dec. J(2000) 


SED Class 




K mag 


L' mag 


04108+2803A 


04 13 53.26 


+28 11 23.6 


Class I 


0.04 


10.32 


9.4 


04108+2803B 


04 13 54.69 


+28 11 33.1 


Class I 


0.52 


10.65 


7.9 


04158+2805 


04 18 57.76 


+28 12 22.7 


Class I 


0.16 


10.66 


9.4 


04181+2654A 


04 21 11.47 


+27 01 09.4 


Flat 


0.24 


10.47 


8.5 


04181+2654B 


04 21 10.39 


+27 01 37.3 


Flat 


0.19 


10.98 


9.4 


04239+2436 


04 26 57.1 


+24 43 36 


Class I 


0.31 


9.88 


7.2 


04260+2642 


04 29 05.03 


+26 49 04.4 


Flat 


0.27 


11.33 


10.3 


04264+2433 


04 29 30.08 


+24 39 55.1 


Class I 


0.18 


11.27 


10.7 


04295+2251 


04 32 32.07 


+22 57 30.3 


Flat 


0.15 


10.05 


8.7 


04361+2547 


04 39 13.87 


+25 53 20.6 


Class I 


0.21 


10.32 


8.4 


04365+2535 


04 39 35.01 


+25 41 45.5 


Class I 


0.22 


10.42 


7.0 



^<7k is the standard deviation of multiple K-band brightness measurements (in magnitudes) of these 
Class I stars, as determined in the study of Park & Kenyon (2002). 
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Table 2. Detected Spectral features 



Species 


Wavelength (/im) 


Note 


Ha S(0) 


1.97 


2 


H2 S(2) 


2.03 


2 


He I 


2.06 


6 


H2 S(l) 


2.12 


2 


HI Brackett 7 


2.16 


1 


Na 


2.20 


3 


Ha S(3) 


2.22 


2 


Ca 


2.26 


3 


CO 


2.29 


3,4 


CO 


2.32 


3,5 


CO 


2.35 


3,4 


CO 


2.38 


3,4 


Ha Q(l) 


2.40 


2 


Ha Q(2) 


2.41 


2 


Ha q(3) 


2.42 


2 


H3O ice 


~3.05//m 


5 


C-H Stretch 


~3.47//m 


5 


HI Pfund 7 


3.74 


1 


HI Brackett a 


4.05 


1 



Note. — 1 - Magnctosphcric emission from 
material accreting onto the star, 2 - Molec- 
ular hydrogen emission, likely from shocks in 
outflows, 3 - Photospheric absorption features 
characteristic to young stars, 4 - CO band-head 
emission from warm inner regions of YSO accre- 
tion disks or winds, 5 - Broad, solid-state water 
ice absorption and shallow C-H stretch absorp- 
tion from ices frozen onto dust grains along the 
line of sight. 6 - Helium I emission from ac- 
cretion driven winds in stars with strong mass 
accretion. 
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Table 3. Spectral Line Fluxes and Detection Limits (xlO W/m^/s) 



Star Name 


Br7 


Bra 


Pf7 


Ha 8(0) 


Ha 8(2) 


H2S(1) 


H28(3) 


H2Q(1) 


H2Q(2) 


H2Q(3) 


04108+2803B 


3.4+0.4 


96+4 


27.7+0.6 


<2.1 (3(t) 


<3 (3<t) 


1.3+0.4 


<2.4 (3(7) 


<3 {3a) 


<3 (3<t) 


<3 {3a) 


04158+2805 


1.8+0.6 


10+2 


<3.6 


<4.5 


<3 


2.7+0.5 


<1.5 


<3 


<3 


<3 


04181+2654A 


8.4+1.2 


53+3 


15+3 


<3 


<3 


3.5+0.9 


<3 


<4.5 


<4.5 


<4.5 


04181+2654B 


<1.8 


<6 


<3 


<4.5 


<3 


<2.4 


<2.4 


<6 


<6 


<6 


04239+2436 


53.2+0.6 


319+4 


117+2 


<3 


<2.4 


5.7+0.6 


<2.1 


<6 


<6 


<6 


04260+2642 


<2.1 


4+1 


<1.5 


8.3+0.9 


3.1+0.9 


39.1+0.5 


<2.1 


11+2 


4+2 


8+2 


04264+2433 


3.5+0.7 


<6 


<4.5 


22+1 


4+1 


9.6+0.5 


4.5+1.5 


12+1 


6+1 


9+1 


04295+2251 


4+1 


18+2 


<3 


6+2 


4+1 


11+1 


5+1 


7+2 


<4.5 


4.5+1.5 


04361+2547 


9.3+0.6 


28.0+0.8 


10+1 


9.5+0.6 


<3 


5.7+0.7 


4+1 


12.6+0.8 


5.5+0.8 


8.0+0.9 


04365+2535 


13+1 


131+2 


27+3 


<3 


<2.1 


2+1 


<2.4 


<3 


<3 


<3 
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Table 4. Water Ice 



Star Name 






Auice 




Nice 


04108+2803B 


640 


3.05 


550+60 


1.49+0.05 


4.0+0.5x101* 


04158+2805 


910 


3.12 


910+90 


0.6+0.1 


2.6+0.8x10^* 


04181+2654A 


750 


3.04 


780+90 


1.26+0.07 


4.9+0.6x10^* 


04181+2654B 


_i 


3.08 


700+100 


2.2+0.2 


8+2x10^* 


04239+2436 


670 


3.05 


510+90 


1.5+0.1 


4+lxlOi* 


04260+2642 


980 


3.06 


630+90 


0.39+0.09 


1.2+0.5x10^* 


04264+24332 


1310 










04295+2251 


890 


3.06 


600+100 


1.0+0.1 


3.0+0.9x10^* 


04361+2547 


790 


3.04 


730+70 


2.2+0.1 


8+lxlOi* 


04365+2535 


560 


3.05 


820+150 


2.43+0.07 


10+2x10^* 



iPor IRAS 04181+2654B, no blackbody temperature could be found that 
had less than ~15% flux deviation from both the K and L-band continuum 
spectra. 

^The L-band spectra of IRAS 04264+2433 was too noisy to identify with 
confidence the optical depth and FWHM of a water ice absorption feature. 
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Table 5. Spectral Synthesis Best- fit Models 



Star Name 


Spectral Type 




TK (2.22 ^jm) 


r^ (3.76 fim) 


Tdust(K) 




04158+2805 
04181+2654A 

04181+2654B 

04260+2642 

04295+2251 


M6±l 
M3tl 
K7+1 
K5+2 
K7±2 


16.3 
17.8 
40.5 
20.7 
28.0 


0.25+0.06 

1.3+0.2 

00+"-"'' 
-().()() 

0.28+0.07 
0.09+0.03 


0.8+0.1 
6.2+0.5 

1.1+0.1 
0.38+0.07 


985 
905 

1020 
1040 


0.33 
1.45 

0.38 
1.12 



Note. — For IRAS 04181+2654B, the models that fit the steep slope of the K-band spectra had 
zero infrared veiling but always overestimate the L-band continuum shape by ~50% (See Figure 
3c). 
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Table 6. Visual Extinction Determined Using Five Independent Methods 



Star Name 


A„ (ice) 


A„ (phot) 


A„(H2) 


(spec) 


A^(Bra:/Br7) 


04108+2803B 


19.1 


16.3 






33.9 


04158+2805 


9.3 


2.7 




16.3 


9.2 


04181+2654A 


16.7 


22.9 




17.8 


12.1 


04181+2654B 


26.6 


40.8 




40.5 




04239+2436 


19.2 


20.0 






10.6 


04260+2642 


7.3 


6.0 


15.7 


20.7 




04264+2433 




9.7 


neg. 




neg 


04295+2251 


14.1 


17.5 


3.3 


28.0 


6.1 


04361+2547 


26.4 


20.8 


9.1 




0.0 


04365+2535 


29.3 


12.0 






18.4 



Note. — "Negative" values of A„ are derived for IRAS 04264+2433 from the 
H2 and Bra/Br7 line ratios, which means that this star is viewed in the infrared 
predominantly in scattered light. 



-34- 



Table 7. Stellar Accretion Luminosity 



Star Name 


Adopted A„ 


Dereddened Br7 Flux 


Lace (L©) 




(Locc/Lsoi) 




(xlO" 


W/m^Z/im) 








04108+2803B 


33.9 


112.0 


0.15 


0.6 


0.25 


04158+2805 


16.3 


9.7 


0.007 


0.2 


0.04 


04181+2654A 


17.8 


52.7 


0.06 


0.3 


0.2 


04181+2654B 


40.5 


<117 


<0.16 


0.3 


<0.53 


04239+2436 


20.0 


386 


0.7 


1.3 


0.54 


01260+2642 


20.7 


<17.8 


<0.01 






04264+2433 


9.7 


9.5 


0.007 


0.3 


0.02 


04295+2251 


28.0 


71.8 


0.09 


0.5 


0.17 


01361+2547 


26.1 


142 


0.20 


2.5 


0.08 


131)5+2535 


29.3 


2(>7 


O.li 


1.9 


0.23 


Additional Data from Muzerolle et al. (1998) 


04361+2547 






0.08 


2.5 


0.03 


04016+2610 






0.03 


3.6 


0.008 


04239+2436 






0.24 


1.3 


0.18 


04263+2426B 






0.05 


5.5 


0.009 


04489+3042 






0.02 


0.3 


0.07 


04292+2422 






0.18 


1.2 


0.15 



Note. — ^ Bolometric luminosities are from Chen et al. (1995). 
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Table 8. Stellar Mass Accretion 



Star Name 1 Myo Mass (Mq) M (Me/yr) log(M) (Mo/yr) 



04158+2805 


0.05 


7.5x10"^ -8.1 


04181+2654A 


0.32 


9.0x10-8 -8.0 


04181+2654B 


0.68 


<2.5xl0-* <-7.6 


04260+2642 


0.89 


<4.5xlO-io <-9.3 


04295+2251 


0.68 


1.2x10-** -7.9 


Additional Data from White 


& Hillenbrand (2004) 


04158+2805 


0.05 


<3.2xl0-i" <-9.5 


04260+2642 


0.89 


1.6x10-'^ -6.79 



Note. — Approximate masses were derived by fitting the temperature 
estimate of the star to the 1 My isocrone using the Siess et al. (2000) 
evolutionary models (White & Hillenbrand 2004). 
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Table 9. Summary of Class I Properties 



Star Name 




ct 






Silicates 


Mf (mq) 


Polarization (Pk) 


04365+2535 


172 


+1.27 


0.23 


2.43 




~0.25 


14.6% 


04108+2803B 


(205) 


(+0.62) 


0.25 


1.49 


Abs+Emis 


0.008U 


1.6% 


04239+2436 


236 


+1.15 


0.54 


1.5 


Abs 


0.09 


7.4% 


04181+2654A 


(346) 


(+0.10) 


0.2 


1.45 


Abs 


(0.12) 


1.9% 


04361+2547 


(144) 


(+1.55) 


0.08 


2.2 




0.25 


21.2% 


04181+2654B 


(346) 


(+0.10) 


<0.55 


2.2 


Emis+Abs 


(0.12) 


1.7% 


04295+2251 


447 


+0.21 


0.21 


1.12 


Emis+Abs 


0.025U 


2.6% 


04158+2805 


528 


+0.72 


0.04 


0.33 




0.025U 




04260+2642 




+0.25 




0.38 




0.025U 


1.8% 


04264+2433 


(252) 


(+0.98) 


0.02 




Emis 


0.007U 


4.8% 



Note. — Values in parenthesis were derived from spatially unresolved fluxes (Chen et al. 1995; White & 
Hillenbrand 2004; Motte & Andre 2001). Accretion luminosities and water-ice optical depths are from this study, 
characteristics of the lO^m silicate profiles are from Kcssler-Siliacci et al. (2005) and Watson et al. (2004). The 
Envelope masses within 4200AU are from Motte & Andre (2001), a U designates sources that were unresolved in 
their millimeter continuum study. Polarization values, Pk , are from Whitney et al. (1997). 
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Figure Captions: 

Figure 1: The near infrared spectra of ten Class I and Transition protostars. The 
upper panel for each star shows the full 2-4/xm range, and the lower inset shows a zoomed 
view of the K-band region (1.96-2.45/xm). Emission and absorption features detected in the 
spectra are identified in the plots and described in Table 2. 

Figure 2: The K-band Spectra of five protostars that have photospheric absorption 
features. Overplotted are best fit models derived by applying a (constant) veiling value 
to a spectral standard star and reddening it using ISM extinction. The model results are 
presented in Table 5. 

Figure 3: The 2-4/im spectra of the 5 Class I protostars that have photospheric absorp- 
tion features. Overplotted are the best fit models derived by adding a wavelength dependent 
infrared excess emission to a spectral template standard, reddening spectrum using an in- 
terstellar extinction law (A(A)=At,(0.55/A)^'^), and finally applying a wavelength dependent 
ice-band absorption model using available laboratory data for pure H2O ice (See §3). The 
unexplained long wavelength wing of the water ice absorption feature is apparent in each of 
the spectra (3.25-3.8/im). For IRAS 04181+2654B (panel c), the best model gave a poor fit 
to slope of the K-band. Overplotted in green is the best fit to the K-band spectra, extended 
to 4.0//m. In panel e) for IRAS 04295-1-2251, overplotted in blue is a fit for a model that has 
a reduced — residual and demonstrates uncertainties in individual model fits. 
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v=l-0 5(1) Line Flux (xlO"^^ W/mVs) 



Fig. 5. — A comparison of the emission lines of H2 Q(3) versus H2 S(l). 

The intrinsic hne ratio is 0.7; deviations from this value are caused by either extinction or 
scattering. Overplotted in the figure are the resulting line ratios from A^=10mag, 20 mag of 
extinction, and a scenario where the molecular hydrogen flux is unobscured but seen purely 
in scattered light with a Rayleigh wavelength dependence. 
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Fig. 6. — A comparison of the emission features of Bra plotted versus Br7. Intrinsic 
line ratios in a Baker & Menzel (1938) Case B treatment are in the range of 2.85 to 3.17 
(Brockelhurst 1971); deviations from this range are presumably caused by either extinction 
or scattering. Overplotted in the figure are the ranges that would result from A^,=10mag 
extinction, 20 mag of extinction, and a scenario where the Brackett fluxes are seen purely in 
scattered hght, as in Figure 5. IRAS 04239+2436 is outside to the upper right of the plotted 
limits. 
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Fig. 7. — The 2-Afim spectra for all 10 targets from this study, plotted from bottom to top 
in order of increasing K-L' color. The lower stars in the plot are best described as transition 
Class I/Class II stars, and the upper are true Class I protostars. Included to the left of the 
figure is a note on the order of the spectra from Figure 1. IRAS 04181+2654B and IRAS 
04295+2251 could be more evolved Class II T Tauri stars that are heavily obscured (see §6). 



